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1. Introduction 

The ground-breaking work of the imag- 
ing atmospheric Cherenkov telescope (IACT) 
Whipple led to the earliest detections of 
sources in the VHE (i.e., -0.1-100 TeV) band, 
the Crab pulsar and nebula (Weekes et al. 
1989) and the blazar Mrk 421 (Punch et 
al. 1992). These pioneering attempts initiated 
VHE astrophysics. Following on, the first gen- 
eration of major IACTs - that included CAT 
(1996-2003), CANGAROO (1992-2001) and 
HEGRA (1993-2002), in addition to Whipple 
(1969-2003) itself - broadened the discovery 
potential of the new field by detecting sev- 
eral more blazars and Galactic sources. Thanks 
to reduced low-energy thresholds, improved 
sensitivities, wider-field cameras, and lighter 
mechanical structures, the current second- 
generation IACTs - i.e., H.E.S.S. (2003-), 
MAGIC (2004-), CANGAROO III (2004-), 
and VERITAS (2006-) - have taken VHE as- 
trophysics into maturity. 

In this paper we highlight some recent 
progress in VHE astrophysics obtained with 
the current generation of instruments. 
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2. Galactic sources 

In a seminal VHE survey of the Galactic 
plane, the southern-located H.E.S.S. telescope 
discovered 14 previously unknown sources 
(Aharonian et al. 2006a). Further Galactic 
sources, accessible from the northern hemi- 
sphere, were subsequently observed with the 
MAGIC telescope (e.g., Albert et al. 2007a). 
Proposed counterparts of such Galactic VHE 
sources include supernova remnants (SNRs), 
pulsar wind nebulae (PWNe), and accreting bi- 
naries. Whatever their detailed nature, it is ex- 
pected that Galactic VHE sources are related 
to evolutionary endproducts of massive, bright, 
short-lived, stellar progenitors. Hence, these 
Galactic VHE sources are immediate tracers of 
the current star formation. 



2.1. Supernova remnants 

Galactic cosmic rays have long been suspected 
to be produced at supernova (SN) shock fronts 
via diffusive acceleration. If the observed VHE 
y-rays were found to be generated through 
the hadronic channel, via 7r° decay follow- 
ing pp interaction with the dense molecular 
clouds embedding the short-lived SN progen- 
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itor, then the acceleration by SNe of nuclei 
to energies of the order of the knee in the 
CR spectrum would be virtually proven (e.g., 
Torres et al. 2003). However, it is difficult 
to disentangle the hadronic VHE component 
from the leptonic one, produced by inverse- 
Compton (IC) scattering of interstellar radi- 
ation field photons (in the inner Galaxy) or 
CMB photons (in the outer Galaxy) off ultra- 
relativistic electrons (e.g., Porter et al. 2006), 
by measuring y-rays over only a decade or 
so in energy. The VHE data of RX J1713.7- 
3946 can be explained in terms of either chan- 
nel, leptonic/hadronic if the relevant magnetic 
field is low/high (fi~10/100^G: Aharonian et 
al. 2006b and Berezhko & Volk 2006). Data 
in the ~0. 1-100 GeV band, such as those 
to be provided by AGILE and GLAST, are 
clearly needed to discriminate between the 
two channels. (For Cas A the high magnetic 
field, B~\ mG, suggests a mostly hadronic 
[Berezhko et al. 2003] VHE emission [Ona- 
Wilhelmi et al. 2007].) 

Whatever the details, the detection of pho- 
tons with energy > 100 TeV from RX J 17 13.7- 
3946 is a proof of the acceleration of pri- 
mary particles in SN shocks to energies well 
above 10 14 eV. The differential VHE spec- 
tral index is ~2.1 all across this SNR, sug- 
gesting that the emitting particles are ubiq- 
uitously strong-shock accelerated, up to en- 
ergies ~200/100 TeV for primary CR pro- 
tons/electrons if the hadronic/leptonic channel 
is at work (Aharonian et al. 2007). This is get- 
ting close to the knee of the CR spectrum, 
~10 15 5 eV, that signals the high-energy end of 
the Galactic CR distribution (e.g., Blasi 2005). 

Circumstantial evidence supports a 
hadronic origin of the VHE emission. In 
several expanding SNRs the X-ray brightness 
profile behind the forward shock is best 
explained as synchrotron emission from 
energetic electrons in high magnetic fields, 
B~O(10 2 )juG, i.e. -100 times larger than typ- 
ical interstellar medium (ISM) values. Such a 
large amplified magnetic field disfavors the IC 
interpretation of the VHE data. Furthermore, 
in the remnant H.E.S.S. J1834-087 the maxi- 
mum of the extended VHE emission correlates 
with a maximum in the density of a nearby 



molecular cloud (Albert et al. 2006a) - which 
suggests hadronic illumination of the target 
molecular cloud. 



2.2. Pulsars and associated nebulae 

Discrimination between different processes of 
pulsar magnetospheric emission (e.g., polar- 
cap vs outer-gap scenario) is one clear goal of 
VHE astrophysics. Polar-cap (e.g., Daugherty 
& Harding 1982, 1996) and outer-gap (e.g., 
Cheng et al. 1986; Romani 1996) models es- 
sentially differ by the location of the gap in 
the pulsar magnetosphere. In the former case 
this is close to the neutron star (NS) surface, 
whereas in the latter it is further away from 
it. Thus, the influence of the magnetic field 
(B~10 n ~ 13 G) is crucially different in these 
models. In polar-cap models, it produces ab- 
sorption (due to y+B — »e ± ) leading to a super- 
exponential cutoff of the emission (mostly cur- 
vature radiation). In outer-gap models, only a 
(milder) exponential cutoff is present, and the 
highest photon energies depend on the elec- 
tron energy: in these models a VHE IC compo- 
nent may thus arise from the upscattering, by 
such energetic electrons, of their emitted syn- 
chrotron photons or of the X-ray photons re- 
leased by NS heating. Both models can deal 
with current observational constraints. 

The current situation is exemplified be the 
cases of the Vela and PSR B 195 1+32 pulsars 
(Aharonian et al. 2006c; Albert et al. 2007b). 
In the latter case upper limits to the pulsed 
emission imply a cutoff energy E c <2>2 GeV. In 
both, Vela and PSR B 195 1+32, IC emission at 
TeV energies as predicted by outer-gap models 
is severely constrained, although not all outer- 
gap models are ruled out. Deeper sensitivities 
can test the models further, and certainly no 
test of these models in the range 10-100 GeV 
has yet been achieved. 

Millisecond pulsars, that have lower mag- 
netic fields (fi~10 8 ~ 10 G), in polar-cap models 
could have E c > 100 GeV (i.e., £>£-'), hence 
their (pulsed) VHE emission could be substan- 
tial. Detecting it would be a test of the polar- 
cap theory (Harding et al. 2005) - however, no 
positive detection has been reported so far. 
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The Crab nebula, a steady emitter that is 
used as a calibration candle, has been observed 
extensively from the radio up to ~70 TeV. 
No pulsed (magnetospheric) VHE emission 
was found in MAGIC data, that implies 
£ e <30 GeV (Albert et al. 2007c). The steady 
nebular spectrum, measured in the ~0.03- 
30 GeV range by EGRET and in the -0.06- 
70 TeV by several IACTs, shows a bump that 
starts to dominate at ~1 GeV and peaks at 
~50 GeV: this component results from IC scat- 
tering, by the synchrotron-emitting electrons, 
of softer photons in the shocked wind region - 
i.e., synchrotron, FIR/mm or CMB photons. In 
spite of its detected IC y-ray emission, how- 
ever, the Crab nebula is not an effective IC 
emitter as a consequence of its high nebular 
magnetic field (5-0.4 mG). 

PWNe - i.e., pulsars displaying a promi- 
nent nebular emission - currently constitute 
the most populated class of identified Galactic 
VHE sources (7 identifications - e.g., Gallant 
2006). The VHE emission of PWNe is likely 
of leptonic origin. Let us consider the case of 
H.E.S.S. J1825-137, for which spectra have 
been measured in spatially separated regions 
(Aharonian et al. 2006d). In these regions, the 
VHE spectra steepen with increasing distance 
from the pulsar, and the VHE morphology is 
similar to the X-ray morphology: furthermore 
the low derived magnetic field (few fiG) im- 
plies that synchrotron X-ray emission is due to 
electrons of energy higher than the y-rays. This 
suggests that the observations can be inter- 
preted in terms of very energetic electrons that 
efficiently lose energy via synchrotron losses, 
aging progressively more rapidly as they are 
farther away from the acceleration site, and 
produce VHE y-ray s via IC scattering. 

2.3. TeV binaries 

In both SNRs and PWNe particle accelera- 
tion proceeds on the parsec distance scales in 
the shocks formed in interactions of either SN 
ejecta or pulsar winds with the ISM. A differ- 
ent population of much more compact parti- 
cle accelerators, which has been revealed by 
current IACTs, is formed by the TeV bina- 
ries (TVBs). These systems contain a com- 



pact object -either a NS or a black hole (BH)- 
that accretes, or interacts with, matter outflow- 
ing from a companion star: hence they are 
VHE-loud X-ray binaries (XRBs). Four TVBs 
have been detected so far: PSR B1259-63 
(Aharonian et al. 2005a), LS 5039 (Aharonian 
et al. 2005b, 2006e), LS I +61 303 (Albert et 
al. 2006b), and Cyg X-l (Albert et al. 2007d). 

PSR B 1259-63 is powered by the rotation 
energy of its young 48 ms pulsar, and strong 
VHE emission is observed from this system in 
pre- and post-periastron phases when the rel- 
ativistic pulsar wind collides with the dense 
equatorial wind blowing from the companion 
Be star. 

LS I +61 303 and LS 5039 may share a 
similar structure. The former is composed of a 
compact object and a Be star in a highly eccen- 
tric orbit. Its VHE emission, whose variability 
constrains the emitting region of LS 1+61 303 
to be larger than the binary system's size, ap- 
pears correlated with the radio emission and 
does not peak at periastron, where M is ex- 
pected to be largest. This picture favors an IC 
origin of VHE emission, as probably the most 
efficient at the relatively large scales of the sys- 
tem at peak emission. 

The emission from Cyg X-l is point-like 
and excludes the nearby radio nebula pow- 
ered by the relativistic jet. Cyg X-l is the first 
stellar-mass BH, and hence the first established 
accreting binary, detected at VHE frequencies. 

2.4. Galactic center 

The possibility of indirect dark matter (DM) 
detection through its annihilation into VHE 
y-rays has aroused interest to observe the 
Galactic Center (GC). H.E.S.S. and MAGIC 
observed the GC, measuring a steady flux con- 
sistent with a differential power-law slope of 
~2.2, up to energies of ~20 TeV with no 
cutoff (Aharonian et al. 2004; Albert et al. 
2006c). Within the error circle of the mea- 
surement of the central source H.E.S.S. J1745- 
290 are three compelling candidates for the 
origin of the VHE emission: the shell-type 
SNR Sgr A East, the newly discovered PWN 
G 359.95-0.04, and the supermassive BH 
Sgr A* itself. Plausible radiation mechanisms 
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include IC scattering of energetic electrons, the 
decay of pions produced in the interactions of 
energetic hadrons with the ISM or dense radia- 
tion fields, and curvature radiation of UHE pro- 
tons close to Sgr A*. These considerations dis- 
favor DM annihilation as the main origin of the 
detected flux, whereas a more conventional as- 
trophysical mechanism is likely to be at work 
(e.g., Aharonian et al. 2006f). Furthermore, 
the lack of flux variability on hour/day/year 
timescales suggests that particle acceleration 
occurs in a steady object, such as a SNR or a 
PWN, and not in the central BH. 

The GC diffuse emission correlates with 
molecular clouds and suggests an enhanced 
CR spectrum in the Galactic center (Aharonian 
et al. 2006g). Its morphology and spectrum 
suggest recent in situ CR acceleration: because 
the photon indexes of the diffuse emission and 
of the central source H.E.S.S. J1745-290 are 
similar, the latter source could be the accelera- 
tor in question. 

3. Star-forming galaxies 

Diffuse y-ray emission from pp interactions of 
CR nuclei with target ISM and photons makes 
up -90% of the >100 MeV luminosity of the 
Milky Way (Strong et al. 2000). However, the 
VHE flux from a galaxy like the Milky Way 
located 1 Mpc away would be well below cur- 
rent IACT sensitivities. Indeed, only loose up- 
per limits on the VHE flux from normal galax- 
ies have been obtained, even for local galaxies 
and for the VHE-bright starburst galaxies (e.g., 
Torres et al. 2004). Detailed models of VHE 
emission from NGC 253 (Volk et al. 1996; 
Paglione et al. 1996; Domingo-Santamana & 
Torres 2005) and for Arp 220 (Torres 2004) 
are only loosely constrained by current upper 
limits (Aharonian et al. 2005c; Albert et al. 
2007e). 

4. Active galactic nuclei 

Supermassive black holes (SMBHs) reside in 
the cores of most galaxies. The fueling of 
SMBHs by infalling matter produces the spec- 
tacular activity observed in active galactic nu- 
clei (AGNs). 



The current AGN paradigm includes a cen- 
tral engine, most likely a SMBH, surrounded 
by an accretion disk and by fast-moving 
clouds, which emit Doppler-broadened lines 
(Un-y & Padovani 1995; Padovani 2007). In 
~10% of all AGNs, the infalling matter turns 
on powerful collimated jets that shoot out from 
the SMBH in opposite directions, likely per- 
pendicular to the disk, at relativistic speeds. 

If a relativistic jet is viewed at small an- 
gle to its axis the observed jet emission is 
amplified by relativistic beaming Q and dom- 
inates the observed emission. Such sources 
are called blazars. Given the blazars' com- 
pactness (as suggested by their short variabil- 
ity timescales), all GeV/TeV photons would 
be absorbed through pair-producing yy colli- 
sions with target X-ray/IR photons. Beaming 
ensures the intrinsic radiation density to be 
much smaller than the observed one, so that y- 
ray photons encounter a much lower yy opac- 
ity and hence manage to leave the source: re- 
versing the argument, y-ray detection is a proof 
of strongly anisotropic (e.g., beamed) emis- 
sion. 

The spectral energy distributions (SEDs) 
of blazars are generally characterized by two 
broad humps, peaking at, respectively, IR/X- 
ray and GeV-TeV frequencies (Ulrich et al. 
1997). Analyses of blazar SEDs (Fossati et al. 
1998; Ghisellini et al. 1998) have suggested 
that: (i) higher/lower-luminosity objects have 
both humps peaking at lower/higher frequen- 
cies (they are called, respectively, LBLs and 
HBLs); (ii) the luminosity ratio between the 
high- and low-frequency bumps increases with 
luminosity; (Hi) at the highest luminosities the 
y-ray output dominates the total luminosity. 



1 Defining the relativistic Doppler factor as 
(5=[r(l-/3cos 8)]~ l (with p=v/c the jet speed nor- 
malized to the speed of light, T= 1 / -yS 2 ), and 
6 the angle w.r.t. the line of sight), the observed and 
intrinsic luminosities at a given frequency / are re- 
lated by Lf s =S p Lf with p~2-3, and the variability 
timescales are related by Af b S =<5 _1 A? cm . For 8~0° 
and 6~2 T the observed luminosity can be amplified 
by factors -400-10 4 (for, typically, r~10 and p~2- 
3); whereas 6~l/T implies 6~T, with a luminosity 
amplification of ~10 2 -10 3 . 
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The mainstream interpretation of the 
blazars SEDs is synchrotron-Compton emis- 
sion, i.e. synchrotron emission (peaked in the 
IR/X-ray range) from a time-varying popula- 
tion of ultra-relativistic electrons moving in a 
strong magnetic field, and IC emission (peaked 
in the ~100MeV-100GeV range) from soft 
photons scattering off energetic electrons. 
Depending on the relative efficiency of the rel- 
ativistic particles' cooling through scattering 
with photon fields that are internal to jet or ex- 
ternal to it, the synchrotron and Compton com- 
ponents peak at, respectively, UV/X-ray and 
GeV-TeV energies (synchrotron-self-Compton 
[SSC] scheme: e.g., Maraschi et al. 1992) or at 
IR/optical and MeV-GeV energies (external- 
IC [EIC] scheme, see Dermer & Schlickeiser 
1993). Hybrid SSC/EIC models have also been 
proposed (Ghisellini 1999). Alternative mod- 
els of VHE emission involve, e.g., two electron 
populations, one -primary- accelerated within 
the jet and the other -secondary- generated 
by electromagnetic cascades initiated by pri- 
mary protons/nuclei that had been accelerated 
in the jet (Mannheim 1993); or a population of 
extremely energetic protons emitting by syn- 
chrotron radiation (Aharonian 2000). 

The emitting particles are accelerated 
within the relativistic jets which transport en- 
ergy from the central SMBH outwards (Rees 
1967). In the popular SSC framework this pro- 
cess is approximated with a series of relativis- 
tically moving homogeneous regions (blobs), 
where particle acceleration and radiation take 
place (e.g., Maraschi et al. 1992). The X-ray 
and y-ray emission, with its extremely fast 
and correlated multi-frequency variability, in- 
dicates that often a single region dominates the 
emission. 

VHE data are of crucial importance to 
close the SSC model. Even in the simplest one- 
zone SSC model of blazar emission, knowl- 
edge of the whole SED up to the VHE regime 
is required for a complete description of the 
emitting electrons' distribution and environ- 
ment (e.g., Tavecchio et al. 1998). The param- 
eters that specify the properties of the emit- 
ting plasma in the basic SSC model are: the 
electron distributions normalization, low- and 
high-energy slopes, and min/break/max en- 



ergy, and the plasma blobs magnetic field, 
size and Lorentz factor. Knowing only the 
IR/X-ray peak would give info on the shape 
(i.e., the broken-power-law slopes) of the elec- 
tron distribution but would leave all other pa- 
rameters unconstrained (e.g., Tavecchio et al. 
1998). However, accurate knowledge of blazar 
emission mechanism(s) requires simultaneous 
broadband y-ray and X-ray (i.e., IC and syn- 
chrotron) data. In fact, a simultaneous SED can 
act as a snapshot of the emitting population of 
particles at a given time. 

Blazar observations have been a top prior- 
ity for VHE astrophysics ever since the discov- 
ery of TeV emission from Mrk 421 (Punch et 
al. 1992). To date, firm blazar TeV detections 
include 15 HBLs and 1 LBL. (One further non- 
blazar AGN, M 87, has also been detected, see 
Aharonian et al. 2006h).) 

The known TeV blazars are variable in flux 
in all wavebands. Even simple one-zone ho- 
mogeneous SSC modeling predicts the X-ray 
and TeV flux variability to be closely corre- 
lated, both emissions being linked to the same 
electron population. Observational evidence, 
although still statistically limited, supports this 
prediction (e.g., Pian et al. 1998). (Some TeV 
flares that show no simultaneous X-ray coun- 
terpart may be explained as IC radiation from 
an additional, very compact and dense, elec- 
tron population - see Krawczynski et al. 2004.) 
Blazar variability, in flux and spectrum, has 
been observed at VHE frequencies down to 
minute timescales. For Mkn 501, observed 
with the MAGIC telescope at > 100 GeV dur- 
ing 24 nights between May and July 2005, the 
integrated flux and the differential photon spec- 
tra could be measured on a night-by-night basis 
(Albert et al. 2007f). During the observational 
campaign, the flux variations (from ~0.4 to ~4 
crab units) were correlated with the spectral 
changes (i.e., harder spectra for higher fluxes), 
and a spectral peak showed up during the most 
active phases. A rapid flare occurred on the 
night of 10 July 2005, showing a doubling time 
as short as ~2 minutes and a delay of ~3 min- 
utes as a function of energy of the emitted pho- 
tons. 

One further aspect of TeV spectra of 
blazars is that they can be used as probes of 
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the Extragalactic Background Light (EBL), i.e. 
the integrated light arising from the evolving 
stellar populations of galaxies (see Hauser & 
Dwek 2001). The TeV photons emitted by a 
blazar interact with the EBL photons and are 
likely absorbed via pair production. Whatever 
its intrinsic shape at emission, after travelling 
through the EBL-filled space, a blazar spec- 
trum will reach the observer distorted by ab- 
sorption. This effect, which is stronger for 
more distant objects (e.g., Stecker et al. 1992), 
is the most likely origin of the avoidance 
zone (i.e., no flat spectra at high redshift) in 
the observed spectral slope vs redshift plot. 
The strength of the absorption is measured by 
Tyy{E), the optical depth for attenuation be- 
tween the blazar, located at a distance D(z), 
and the Earth (Fazio & Stecker 1970; Stecker 
et al. 1992). Usually, either (i) the shape and 
intensity of EBL(z) is assumed, and the TeV 
spectrum is corrected before the SSC mod- 
eling is performed (e.g., de Jager & Stecker 
2002; Kneiske et al. 2004); or (ii) based on 
assumptions on the intrinsic VHE spectrum, 
EBL(z) is solved for: e.g., based on analysis 
of the observed hard VHE spectra of the dis- 
tant blazars 1ES 1101-232 and H 2359-309, 
a low EBL energy density at z~ 0.2 has been 
derived (Aharonian et al. 2006i). The two ap- 
proaches can be used in combination to esti- 
mate the distance to the VHE source (Mazin & 
Goebel 2007). 

5. Gamma-ray bursts 

There is a prevailing consensus that the ba- 
sic mechanism of GRB emission is an expand- 
ing relativistic fireball (Rees & Meszaros 1992, 
Meszaros & Rees 1993, Sari et al. 1998), with 
the beamed radiation (<5~O(10 2 ) due to inter- 
nal/external shocks (prompt/afterglow phase, 
respectively). If so, the emitting particles (elec- 
trons and/or protons) are accelerated to very 
high energies. 

In the fireball shock framework, several 
models have predicted VHE emission during 
both the prompt and afterglow phases of the 
GRB (e.g., Meszaros 2006). This can occur as 
a result of electron self-IC emission from the 
internal shock or the external forward/reverse 



shock. Seed photons can be produced locally 
(through synchrotron, or be the leftover of the 
initial radiation content responsible for the ac- 
celeration of the fireball) or can be produced 
in, e.g., the shell of a previously exploded SN. 
In the latter case, the SN photons may also 
act as targets for the yy absorption, and in 
this case the VHE emission could be severely 
dimmed. If the emission processes are indeed 
synchrotron and IC, then a blazar-like SED is 
predicted, with a double-peak shape extending 
into the VHE band. In such theoretical free- 
dom, VHE observations of GRBs could help 
constraining GRB models. 

MAGIC observed part of the prompt- 
emission phase of GRB 050713a as a response 
to an alert by the Swift satellite (Albert et 
al. 2006d) However, no excess at >175 GeV 
was detected, neither during the prompt emis- 
sion phase nor later - but the upper limits 
to the MAGIC flux are compatible with sim- 
ple extrapolations of the Swift r^l.6 power- 
law spectrum to hundreds of GeV. In general, 
however, the cosmological distances of these 
sources prevent VHE detection (see Albert 
et al. 2007g): the average redshift of the 
GRBs for which MAGIC was alerted (and 
whose z are known) is z=3.22, whereas at 
70 GeV the cosmological y-ray horizon is 
Z~l. Complementary Whipple data (Horan et 
al. 2007) and MILAGRO data (Abdo et al. 
2007) provide upper limits on, respectively, 
the late VHE emission (~4 hr after the burst) 
from several long-duration GRBs, and on the 
prompt/delayed emission from several, reput- 
edly nearby (z~ 0.5), short-duration GRBs. 

6. Dark matter 

Evidence for departure of cosmological mo- 
tions from the predictions of Newtonian dy- 
namics based on visible matter, interpreted as 
due to the the presence of DM, are well estab- 
lished - from galaxy scales (e.g., van Albada et 
al. 1985) to galaxy-cluster scales (e.g., Sarazin 
1986) to cosmological scales (e.g., Spergel et 
al. 2003). 

DM particle candidates should be weakly 
interacting with ordinary matter (and hence 
neutral). The theoretically favored ones, which 
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are heavier than the proton, are dubbed 
weakly interacting massive particles (WIMPs). 
WIMPs should be long-lived enough to have 
survived from their decoupling from radiation 
in the early universe into the present epoch. 
Except for the neutrino, which is the only DM 
particle known to exist within the Standard 
Model of elementary particles (with a relic 
background number density of ~50 ctrT 3 
for each active neutrino species) but which 
is too light (m v ~ 1 eV) to contribute signifi- 
cantly to Q m given the current cosmological 
model, WIMP candidates have been proposed 
only within theoretical frameworks mainly mo- 
tivated by extensions of the Standard Model 
of particle physics (e.g., the R-parity conserv- 
ing supersymmetry [SUSY]). Among current 
WIMP candidates (see Bertone et al. 2005), the 
neutralino, which is the lightest SUSY particle, 
is the most popular candidate. Its relic density 
is compatible with WMAP bounds (see Munoz 
2004). 

WIMPs could be detected directly, via elas- 
tic scattering vith targets on Earth, or indi- 
rectly, by their self-annihilation products in 
high-density DM environments. DM annihila- 
tion can generate y-rays through several pro- 
cesses. Most distinctive are those that re- 
sult in mono-energetic spectral lines, xx~ ^77' 
XX~ *jZ or XX~*yh. However, in most models 
the processes only take place through loop di- 
agrams; hence the cross sections for such final 
states are quite suppressed, and the lines are 
weak and experimentally challenging to ob- 
serve. A continuum y-ray spectrum can also 
be produced through the fragmentation and 
cascades of most other annihilation products. 
The resulting spectral shape depends on the 
dominant annihilation modes (see Bergstrom 
& Hooper 2006), whereas the normalization 
depends on the WIMP's mass and velocity- 
averaged annihilation cross section as well as 
on the DM density profile. 

Once the astroparticle model has been 
chosen (e.g., Bergstrom et al. 1998), the 
biggest uncertainties are of astrophysical na- 
ture. Superposed to any VHE emission from 
the decaying DM (cosmological, non-baryonic 
signal), galaxies can display a VHE emission 
from astrophysical sources associated with the 



visible matter distribution (astrophysical, bary- 
onic signal). The ratio of the former to the latter 
is maximized in small, low-L, low-SFR galaxy. 
This is because the dark-to-visible mass ra- 
tio as well as the central DM density in- 
crease with decreasing luminosity (Persic et 
al. 1996). Clearly, distance dilution of the sig- 
nal opposes detection, so galaxies candidate 
for indirect DM detection should be chosen 
among nearby objects. In conclusion, the best 
obserational targets for DM detection are the 
Milky Way dwarf spheroidal galaxies (e.g., 
Draco, Sculptor, Fornax, Carina, Sextans, Ursa 
Minor). A further issue, stemming from the p 2 
dependence (as a result of annihilation) of the 
normalization integral of the y-ray emission, 
concerns the shape of the inner halo profile, 
i.e. whether the latter is cuspy or cored. Cuspy 
profiles are produced in cosmological N-body 
simulations of halo formation (Navarro et al. 
1997), whereas cored profiles are suggested 
by the measured rotation curves of disk galax- 
ies (Borriello & Salucci 2001) - also in low- 
surface-brightness galaxies, where the local 
self-gravity of baryons is virtually negligible 
(deBloketal.2001). 

These considerations (and uncertainties) 
have been incorporated in detailed predictions 
of the y-ray flux expected from the annihila- 
tion of the neutralino pairs. Outlooks for VHE 
neutralino detection in Draco by current IACTs 
are not very promising: for a neutralino mass 
m^=100 GeV and a variety of annihilation 
modes, and in the favorable case of a maxi- 
mal (cuspy) inner halo profile, VHE detection 
(by MAGIC in 40hr observation) can occur if 
average value of the neutralino's cross section 
times velocity is <<xv>> 10~ 25 cm 3 s _1 , which 
is somewhat larger than the maximum value 
for a thermal relic with a density equal to the 
measured (cold) DM density (but may be fine 
for non-thermally generated relics) in the al- 
lowed SUSY parameter space. The prospects 
are better in the HE range (100 MeV-10 GeV): 
for a maximal (cuspy) halo, 1 yr of GLAST 
observation should be able to yield a detection 
if m x < 500 GeV and <crv>~3xl0" 25 cmV 1 
(Bergstrom & Hooper 2006). 

No evidence of DM annihilation y-rays has 
been unambiguously claimed so far. An ap- 
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parently extended signal from the direction of 
NGC 253 has been definitely interpreted as due 
to hardware malfunction (Itoh et al. 2007). 
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